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Abstract: We have in earlier papers presented an extension of Mueller's dipole cascade 
model, which includes subleading effects from energy conservation and running coupling 
as well as colour suppressed effects from pomeron loops via a "dipole swing" . The model 
was applied to describe the total cross sections in pp and 7*p collisions. 
In this paper we present a number of improvements of the model, in particular related to 
the confinement mechanism. A consistent treatment of dipole evolution and dipole-dipole 
interactions is achieved by replacing the infinite range Coulomb potential by a screened 
potential, which further improves the frame-independence of the model. 
We then apply the model to elastic scattering and diffractive excitation, where we specif- 
ically study the effects of different sources for fluctuations. In our formalism we can take 
into account contributions from all different sources, from the dipole cascade evolution, the 
dipole-dipole scattering, from the impact-parameter dependence, and from the initial pho- 
ton and proton wavefunctions. Good agreement is obtained with data from the Tevatron 
and from HERA, and we also present some predictions for the LHC. 
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1. Introduction 

In high energy pp scattering the cross section for parton-parton subcollisions becomes 
larger than the total cross section. This means that on average there are more than one 
subcollision in a single event, and it was early suggested that hard subcollisions dominate 
the features of high energy scattering and are the cause of the rising cross section. This is 
also the basic assumption in the model by Sjostrand and van Zijl [1] implemented in the 
Pythia event generator, which is able to describe many features of high energy collisions. 

That perturbative dynamics dominate high energy collisions is also supported by the 
large intercept of the BFKL pomeron. Via unitarity and the AGK cutting rules, the 
large subcollision cross section and high probability for multiple collisions have also strong 
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implications for diffraction. Hard diffraction was first observed by the UA8 collaboration 
at the CERN SppS collider [2], and has later been studied in much more detail both at the 
Tevatron (see e.g. [3,4]) and at HERA (see e.g. [5,6]). 

Multiple collisions, unitarity and saturation are conveniently studied in terms of 
dipoles in transverse coordinate space. The dipole model by Golec-Biernat and WusthofF 
(GBW) [7,8] has successfully described both F2 and diffraction in DIS. Mueller's dipole 
cascade model [9-11] reproduces the leading log (linear) BFKL equation, and includes also 
multiple collisions and satisfies the unitarity constraint. The multiple collisions correspond 
to pomeron loops. Mueller's model includes, however, only such loops which axe cut in the 
particular Lorentz frame used in the calculation, but not loops which are fully contained 
in one of the individual dipole cascades. Many attempts (see for example [12, 13] and 
references therein) have been presented including e.g. 2 ^ 1, 2 ^ 4 or more complicated 
dipole vertices, but so far no explicitly frame-independent formalism has been presented. 

As discussed in some detail below, an important part of the NLL corrections to the 
BFKL equation are related to energy conservation. In a series of papers [14, 15] we have 
developed Mueller's model to include both effects of energy-momentum conservation and 
effects of pomeron loops and saturation inside the cascade evolution via a 2 ^ 2 dipole 
transition, called a dipole swing. The swing does not reduce the number of dipoles, rather 
the saturation effect is achieved as the "new" dipoles arc smaller, and therefore have smaller 
cross sections. Although not explicitly frame independent, the numerical result is almost 
independent of the frame used for the calculations. With a simple model for the proton 
in terms of three dipoles, the Monte Carlo implementation also reproduces the total cross 
section both for DIS at HERA and for pp scattering from ISR energies to the Tevatron. 

In this paper we will first make some technical improvements related to confinement, 
and then use the model to study diffractive scattering at HERA and the Tevatron. The 
perturbative calculation has some problems in the IR region, especially with a running 
coupling, and a cutoff for large dipoles is essential for the frame independence and the good 
agreement with data. We here propose to treat this effect of confinement by everywhere 
replacing the Coulomb colour-electric potential by a screened Yukawa potential. 

Both the screening length and the size of the initial proton wave function are deter- 
mined by the confinement mechanism. In the MC implementation these two quantities are 
assumed to be of the same size. This implies that for pp-coUisions and DIS at high the 
model has only two tunable parameters; besides Aqcd only the confinement scale denoted 
by Vmax- For smaller the result is, however, also sensitive to the quark masses in the 
virtual photon wave function. The value of rmax turns out to be very important in order 
to obtain the correct normalization for atot in PP collisions, but the increase of atot with 
the center of mass energy is found to be much less sensitive to this parameter. Once rmax 
is fixed to obtain the correct normalization for atot in PP scattering, the DIS cross section 
is obtained without any further changes. 

Our treatment of elastic scattering and diffractive excitation is based on the cikonal 
approximation and the Good and Walker picture [16]. The result is determined by the 
fluctuations in the collision process originating from the initial wave functions of the proton 
and the virtual photon, from the dipole cascades and from the dipole-dipole scattering 
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probability. In our formalism all these different components give important contributions. 
One result of this is that the impact parameter profile is less steeply falling, i.e. less black 
and white and more "grey", than in models where the dominant fluctuations are assumed 
to come from fluctuations in the impact parameter, b. 

The distribution in the mass, Mx, of the diffractive state can be obtained by a study of 
the collision in different Lorentz frames, as discussed by Hatta et al. [17]. (It is here essential 
that we have a frame-independent formalism.) However, in addition to the fluctuations 
included in this reference and in the GBW approach, we also include fluctuations in the 
evolution of the proton target. 

In section ^ we review briefly the dipole cascade model, discuss the modification of the 
confinement effect, and demonstrate the frame independence of the model. In section ^ we 
discuss the formalism for elastic and diffractive scattering, and the effects of the different 
sources of fluctuations in the collision process. Our results are presented in section |^, and 
the conclusions in section |5[ 

2. Dipole Model and Frame Independence 

We will in this section briefly discuss the cascade model presented in refs. [14, 15], de- 
scribe the modified treatment of confinement, and demonstrate the frame independence by 
showing some quantitative examples. 

2.1 Mueller's cascade model 

The model is based on Mueller's dipole formalism [9-11] in which the small-x evolution is 
interpreted in terms of a dipole cascade. The probability per unit rapidity Y that a dipole 
{x,y) emits a gluon at transverse position z is here given by 



and the evolution of the cascade agrees with the leading order BFKL evolution. As a 
consequence, the total number of dipoles grows exponentially. This also implies a strong 
growth for the total cross section which, however, is tamed by taking multiple dipole 
interactions into account. The scattering probability between two elementary colour dipoles 



Since Mueller's model is formulated in transverse coordinate space, multiple scatterings 
can be included in an eikonal approximation, and a unitarised expression for the total 
scattering amplitude can be obtained as 





with coordinates 







(2.3) 
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2.2 Energy conservation 

The fast growth in leading order BFKL is much reduced by NLL effects. As discussed in 
ref. [18], besides effects from the running couphng (which are also included in our simula- 
tions) the NLL corrections to BFKL evolution are related to the non-singular terms in the 
g gg splitting function, and to the so called energy scale terms, and both these effects 
are related to the conservation of energy and momentum. It should be noted, however, that 
although these NLL corrections to the evolution kernel are essential for the conservation 
of energy and momentum, they are not sufficient to guarantee global energy-momentum 
conservation, which is most easily achieved in a computer simulation [19]. 

The LL BFKL result is obtained from the singular term ~ 1/z in the gluon splitting 
function, which dominates for small z. The recoil for the gluon taking the fraction (1 — z) 
is here neglected. However, for large z-values energy conservation implies that this recoil 
cannot be neglected. (More precisely it is conservation of the positive light-cone component, 
as in a boost invariant formalism the variables XBj and z are defined as fractions of A;+ = 
ko + ki^.) If this suppression can be described by an effective cutoff z < a (with a = 0(0.5)), 
an integral over z is replaced as [20] 



As discussed in ref. [21] the effect of the non-singular terms in the splitting function is just 
such a suppression of large z- values, with Inl/a = 11/12. 

The so called energy scale terms are needed in order for a non-Zc^-ordered ladder to 
look the same if evolved from the bottom up or from the top down. Thus a step down in k± 
must look exactly like a step up in fc^, if generated in the opposite direction. This implies 
that the cascade must be ordered not only in A;+ (via the constraint z < 1) but also in 
the negative light-cone momentum k^. This constraint is therefore also related to energy- 
momentum conservation, and it corresponds to what is commonly called the consistency 
constraint. As discussed in refs. [20,22], it has the consequence that the pole 1/(1 — 7) in 
the characteristic function x(7) is replaced by 1/(1 -|- a; — 7). The minimum of x{l) is then 
shifted to 7 = (1 +uj)/2, and thus a distance a;/2 further away from the poles. Expanding 
this effect to lowest order in a gives the energy scale terms presented in ref. [18]. 

Using the similarities between Mueller's cascade model and the Linked Dipole Chain 
(LDC) model [20,22], we implemented energy-momentum conservation in the dipole cas- 
cade in ref. [14]. The p± of the partons was here associated with the dipole sizes in 
coordinate space. Although the number of dipoles still increases exponentially, the growth 
is significantly reduced and the onset of saturation is delayed. 

2.3 Initial proton and photon wave functions 
2.3.1 Photon 

The splitting 7* qq can be calculated perturbatively, and we use the well known leading 




(2.4) 
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order results for longitudinally and transversely polarized photons: 



r, Q2)|2 = ^ e2 I [^2 ^ _ ^)2](^(^ _ ^)q2 ^ ^2)^2 (^^^^(1 _ ^)Q2 + ^2 

Ko^ (r-^z(l - z)Q2 + | . (2.5) 

Here z is the negative light-cone momentum fraction of the photon carried by q, and r is 
the transverse separation between q and q, and we consider four active quark flavours, with 
an effective light quark mass of 60 MeV and a charm mass of 1.4 GeV, as described in [23]. 



Although the wave functions in eq. (2.5) are well known, there are still ambiguities in 



the initial phase of the evolution of the dipole cascade. This problem is an unavoidable 
consequence of the difficulty to reconcile the fundamentally quantum mechanical process 
with the semiclassical approximation represented by the cascade evolution. 

For a given W in DIS, the total rapidity range available for final state particles is given 
by y = InW'^ /ttiq, where uiq is of the order of^ 1 GeV. We use a Lorentz frame such that 
an interval Yq is on the photon side and the remaining interval y — Yq is on the proton 
side. The kinematics is illustrated in figure ||, which shows the phase-space diagram for 
a DIS event in the (Inp^, y)-plane. Here the positions of the q and q are also indicated. 
Their distance in rapidity from the photon end are given by \np±/z and lnp_|_/(l — z) 
respectively, where p± is the transverse momentum of the quark and the antiquark. 

In the simulations we have assumed that gluon emission is possible only at rapidities 
larger than the rapidity of both q and q. We also identify p± with its typical value 2/r. 
As can be seen from figure |l|, this implies that the interval allowed for the photon-initiated 
cascade is given by 

Y^^ff = Yo- Inf — .) =Yo- -) . (2.6) 

' \mm(z, 1 — z) J \rmm[z, 1 — z) J 

Thus Y^-^'^ depends on both variables r and z. The total range for the evolution 

Yp + Y:;ff = \nW^- Inf . (2.7) 

' \rmm(z, \ — z) J 

may be larger or smaller than Inl/x ~ Inl^/'^ — \nQ^. One could also imagine other 
choices, but we want to emphasize that the difference is sub leading in InM^^, and therefore 
the optimal choice cannot be determined from QCD with present techniques. 

Anticipating the discussion in the following sections, we want in connection to figure 
|l] mention that a calculation of the diffractive cross section corresponds to events which 
have a gap at y = in this particular Lorentz frame. This means that the diffractively 
excited photon is confined within the rapidity range Yq corresponding to a mass limited by 



^In the following the scale mo = 1 GeV will be omitted in writing logarithms. 
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Figure 1: The available phase space for gluon emission in DIS. We have here for simplicity 
assumed that z > 1 — z. The vertical dashed line labeled y = denotes the Lorentz frame in which 
the collision is studied. 

2.3.2 Proton 

The initial proton wavefunction can obviously not be determined by perturbation theory, 
but has to be specified by a model. Our assumption is that the initial proton can be de- 
scribed by three dipoles in a triangular configuration, where the corners could be associated 
with the three valence quarks. Wc admit that this model has deficiencies, and it should 
not be used to study the particle distribution in the proton fragmentation regions. The 
corner connecting two dipoles in a triangle corresponds to an octet colour charge, and not 
a quark triplet, and an alternative model could be three dipoles connected by a junction 
in a Y-shaped configuration. This would imply the need for extra assumptions about the 
behaviour of the colour-neutral junction, and we also believe that the enhanced radiation 
from the three colour octets might mimic radiation from gluon contributions present in the 
proton structure also at low virtuality. 

We have tried two different triangular shapes. In the first one the triangle is assumed 
to be equilateral. In this case the size of the dipoles is assumed to be distributed according 
to a Gaussian 



We found in rcf. [15] that this simple model gives a very good agreement with data on cross 
sections for both pp scattering and DIS. To test the sensitivity to this simple assumption, 
we have now also studied a model where the triangle is not equilateral, but has an arbitrary 




(2.8) 
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10^ 10^ 10"^ 

Vs (GeV) 



Figure 2: The total pp cross section as a function of collision energy for various values of r, 
and Aqcd in units of GeV^^ and GeV respectively. 



shape. Also here the distribution is given by a Gaussian in the sizes of the three sides, ri, 
r2 and r^, 

— - — 2^ -]S{ri + r2 + r3). (2.9) 

The results of the two models are very similar, and we will therefore in the following mainly 
present results obtained with the simple model given by dPi{r). 

The parameter r^ax, which determines the initial dipole size in the proton, is here 
assumed to be the same as the confinement scale in the cascade evolution and the dipole- 
dipole scattering (see section |2]^ below). Along with Aqcd, it is one of the essentially two 
free parameters of our model. We note, however, that the variation of these two parameters 
have similar effects. Thus an increased value for rmax can be compensated by a reduced 
value for Aqcd (and vice versa), leaving the cross sections unchanged as seen in figure |2|. 
We note in particular that the energy dependence is rather insensitive to the parameters 
chosen. 

The probability for a three-body system to contract to a single point should be zero, 
and we have therefore also tested a wavefunction where small r-values are suppressed but 
where (r^) has the same value. We then find essentially the same results for total cross 
sections, but the reduced fluctuations imply that the cross sections for elastic scattering 
and single diffractive excitation become larger. This feature will be further discussed in 



sections 3.2.2 and 



2.4 Dipole swing 

One problem with Mueller's model is the fact that saturation effects are not included inside 
the individual dipole cascades. Thus only pomeron loops which are cut in the particular 
Lorentz frame used are taken into account, and the result is not frame independent. Non- 
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linearities due to multiple interactions are included, but the evolution of the dipole cascades 
obey the linear BFKL equation. 

We argued in ref. [15] that the missing saturation effects can be taken into account 
by including the so called dipole swing in the evolution. Although there is no analytical 
proof that this would give a frame independent formalism, numerical simulations in a MC 
implementation showed that the resulting evolution is approximately frame independent. 



We will show this in more detail in subsection |2.6 . 

The swing is a process in which two dipoles {xi,yj) and {xj,yj) are replaced by two 
new dipoles {xi,yj) and {xj^yj). The process can be interpreted in two ways. The quark 
at Xi and the antiquark at y^ form a colour singlet with probability In this case 

the best approximation of the quadrupole field ought to be obtained by the closest charge- 
anti-charge combinations. Here the swing is therefore naturally suppressed by N^, and it 
should be more likely to replace two given dipoles with two smaller ones. Secondly, we may 
see it as the result of a gluon exchange between the dipoles, which results in a change in 
the colour flow. In this case the swing would be proportional to , which again is formally 
suppressed by A'^^, compared to the splitting process in eq. ( |2.1|) , which is proportional to 

In the MC implementation each dipole is randomly given one of N"^ possible colour 
indices. Only dipoles with the same colour index can swing, and the weight for a swing 
{xi^yi)^ {x2,y2) (3^1,^/2)) (^2,2/1) is determined by factor proportional to 

{xi-yif{x2-y2)'^ 



{xi - y2Y{x2 - yi 
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(2.10) 



This implies that the swing favours the formation of smaller dipoles. The number of dipoles 
is not reduced by the swing, but the fact that smaller dipoles have smaller cross sections 
gives the desired suppression of the total cross section. 

The swing is instantaneous in Y for both the colour multipole and gluon exchange 
interpretations. It is therefore not a vertex in the sense of the dipole splitting whose 
probability is proportional to AY. In the MC implementation the swing is formulated as if 
it was proportional to AF, but its strength is adjusted so that it is effectively instantaneous. 
In ref. [24] it is shown that combining the dipole splitting and the dipole swing, one can 
reproduce all colour correlations induced from the multiple dipole interactions. In case 
each dipole is restricted to single scattering, one can combine a splitting with one swing at 
a time to reproduce all correlations, but without this restriction, the maximum number of 
simultaneous swings needed in combination with a splitting, for a system consisting of N 
dipoles, is — 1. 



2.5 Consistent treatment of confinement 

As the dipole model is formulated within perturbative QCD, confinement effects are natu- 
rally not included. Obviously, one cannot let the dipoles in the cascade become too large, 
and it is then natural to introduce a scale, such as our r^ax parameter, so that large dipoles 
are suppressed. 
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Similarly, the scattering of dipoles is calculated perturbatively, and the interaction 
range is therefore longer than what we would expect from confinement. The formula for 
fij in eq. ( |2.2[ ) is just the two dimensional Coulomb potential which for large distances 
behaves as 

fc~f '"--''-';!r' (2.11) 

where b = ^{{xi + yj — {xj + yj)) is the impact parameter of the dipole-dipole collision. 
Thus the scattering probability falls off only as a power of b, and not as an exponential as 
one would expect from a confining field. The expression for fij can be written as 

f{xi,yi\xj,yj) = ^{A{xi - xj) - A{xi - y^) - A(y^ - x^) + A(y, - yj)f (2.12) 

where A(r) is the Green's function given by 

/■ J2r. ik-r 

To take confinement into account we can replace the infinite range Coulomb potential 
with a screened Yukawa potential. This implies that the Coulomb propagator in 
eq. ( p. 131) is replaced by + M^), where M = 1/r 

max is the confinement scale. In this 

case the integral in eq. ( |2.13| ) is replaced by 



d^k e'^-f 1 



KoirM) (2.14) 



(2^)2fc2 + M2 2tt' 

where Kq is a modified Bessel function. The expression in eq. (2.2) is then replaced by 

\^ 
2 



fij TT ( ^o{\Xi - Xj\/rmax) - Ko{\Xi - yj\/rmaxj - 

2 



-Ko{\yi- Xj\/rmax) + Ko{\yi- yjl/rmax)^ ■ (2.15) 

For small separations where r « Vmax, the function Ko{r /rmax) behaves like ln{rmax/r) 
and we then immediately get the result in eq. ([2.21) . For large separations where r » r^nax, 
Koir /r^ax) fells off exponentially ~ ^^j^Er^E^g-r/rmax ^ expected from confinement. 

To be consistent we should then also modify the dipole splitting kernel accordingly. 



The dipole splitting probability in eq. (2.1) can be written in the form 



dY ~ 271^^ {x - z)^{z - y)^ ~ 2n \{x - zy {y - zf ' ' ^ 



The two terms in this expression are each obtained from the integration 

d^k ke'^ "" _ r d^k e'^ '^ 



I 



{2n)H k^ J (27r)2 k 



(2.17) 
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Figure 3: The quantity Acr/cr = {(7{Yo,Y) - a{Y /2,Y)) / a{Y /2,Y) plotted as a function of 
5 ~ YqIY for different collision energies: full line ^/s = 0.4 TeV, dashed line %fs — 1.8 TeV and 
dotted line -y/s = 3.0 TeV. Also shown is the result excluding the dipole swing for ^fs — 1.8 TeV 
(dash-dotted line). 

Once again making the change — > + M^), and noting that ^Ko{r/rmax) = 

- r.Z,^^ ^iir/rmax), we may replace eq. ( p.l6|) by 

^ -^d^zf—^ ^ _^ Ki{\x - z\/rmax) — 1^ _ ^, Ki{\y - z\/rmax)] ■ 

(2.18) 

For small arguments Ki{r /vmax) ~ ''"j:""' , from which we get back the result in eq. (|2.1(]| ), 
while for large arguments Ki(r /rmax) ~ ^J^EZ^Z^-'^h'm.ax ^ once again we obtain an 
exponentially decaying field. 

2.6 Frame independence 

We will in this subsection demonstrate the frame independence by showing some explicit 
results obtained using the MC implementation. 

The cross section obtained when the right-moving (left-moving) cascade is evolved a 
rapidity distance Yq — ^o) is denoted g^Yq^Y)^ and in figure ^ we show the relative 
difference Aa/a = {ai^Q^Y) - a{Y /2,Y)) / (t{Y /2,Y) plotted vs. 6 = Yq/Y. The figure 
shows results both including the dipole swing and without the swing. Without the swing 
the cross section is too large when (5 — > or (5 — > 1. As expected, the degree of frame 
dependence is increasing for larger t/s, when the saturation effects within cascades become 
more important. When we include the swing, we see that the cross section is (within errors) 
independent of the Lorentz frame used. 

Figure |^ shows similar results for DIS. Here the cross section is not exactly frame 
independent, with a tendency of getting larger as we give more of the total rapidity interval 
to the evolution of the photon. (This may not be very clear from the figures shown but it 
can be seen more clearly for dipole-proton scattering where we do not have ambiguities in 



0.4 TeV 
1 .8 TeV 
3.0 TeV 
no swing 1 .8 TeV 
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Figure 4: The quantity Aa/a = {a{Yo,Y)) - (7{0MY, Y)/a{OMY, Y), plotted as a function of 
S = Yq/Y in DIS for W = 220 GeV and = 8 GcV^ (solid line) and = 14 GeV^ (dashed line). 
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5 

Figure 5: The quantity Aa/a (as defined in figure |^) as a function of 6 for different values of 
collision energy, rmax and Aqcd (in units of TeV, GeV""'^ and GeV respectively) in pp scattering. 

choosing the Y interval as in DIS.) There seems to be two causes for this behaviour, and 
neither of them are related to saturation or the dipole swing. The first cause is the running 
of the coupling, especially when the dipole-dipole scattering amplitudes are calculated 
(see next subsection). The second cause is the treatment of energy conservation in the 
dipole-dipole scattering. In case we use a fixed coupling when calculating the scattering 
amplitude (but still using a running coupling in the evolution), the result appears to be 
more frame independent. These issues need to be handled more carefully, and we intend to 
look at them in forthcoming publications. Note however that the difference is less than 10 
percent for both values, and it is also not dependent on the energy. As our model is not 
expected to have a better accuracy than this, we can consider the model to be essentially 
frame independent also for DIS. 





- 11 - 



< 




Figure 6: The quantity /S.a/a (as defined in figure 0) as a function of 5 in DIS for W — 220 GeV 



for different values of 



and Aqcd (in units of GeV^, GeV ^ and GeV respectively). 



In figures |5| and ^ we also show the same diagrams for different values of the parameters 
rmax and Aqcd, and we can see a similar behaviour as in the previous figures. The 
frame independence of the model is therefore not dependent on the precise values of these 
parameters. 

Naturally also the elastic cross section must be frame independent. Diffraction and 
elastic scattering will be studied in sections |3| and and the numerical result for aei is 
indeed independent of the Lorentz frame used. 



2.7 Running Coupling 

In our simulations we use a running Os both in the dipole splitting and in the dipole-dipole 
scattering probability. Recent NLO studies of the dipole evolution [25,26] have revealed 
a fairly complicated structure for the running of ctg- However, in [27] this was shown 
to simplify in the strongly ordered limits, implying that the relevant scale in the dipole 
splitting is determined by min(r, ri, r2), where r is the mother dipole which splits into ri 
and r2- This is also the scale we have been using in our simulations. 

For the dipole-dipole scattering, the situation is more complicated. We have two 
powers of and there are six different scales involved. With the two colliding dipoles 



and (xj^Uj) we have besides, their sizes, also the distances \xi — Xj\, \xi — 



and l^j — yj\ (cf. eq. ( |2.12| )). Here we have tried two alternatives. In the first case 



the scale of Og for both powers is determined by min(|a;j 
In the second case, we associate one as to each of the colliding dipoles and choose the 



' 2/i I ' I •'^i ■ 



yj\Axi-yj\,\yi 



yi 



yjl^Wi-^jlAVi-yjl) and mm{\xj-yj\,\^^-^j\,\ 
yj\) respectively. To avoid divergencies, as is frozen below the scale 



\yi 



scales min(|a; 

yjlAyi ~ I 

P± = 2/rmax- 

It turns out that the degree of frame dependence in 7*p is similar in both cases and 
the results presented in this paper have all been calculated using the first alternative. 
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3. Diffractive and elastic scattering in the dipole model 
3.1 Formalism 

In this section we will describe the formulas we are going to use in calculating the various 
diffractive and elastic cross sections. We shall rely on the dipole version of the Good and 
Walker picture of diffraction [16] where the scattering eigenstates are given by the dipole 
states. The identification of the QCD parton states as the eigenstates of diffraction is due 
to the work of Miettinen and Pumplin [28] . The situation is complicated by the fact that 
the states of the proton or the virtual photon depend on the Lorentz frame used, and we 
will here quite closely follow the formalism presented in [17]. 

In the Good and Walker picture of diffraction there is a normalized and complete set of 
real particle states {|iV)} with fixed quantum numbers. In addition we have eigenstates of 
the scattering, {|n)}, which also form a complete set of normalized states. Assume that we 
have two incoming particles, one right-moving particle |-R) and one left-moving particle \L). 
These particles can then be diffracted onto the various particle states {|iV)} and {|M)}, 
which carry the quantum numbers of \R) and \L) respectively. The incoming wave is given 

by 

\^i) = \R,L) = Y,c^c^\n,m). (3.1) 

The scattered wave is obtained by operating with ImT on where T is the scattering 
operator. It reads 

\iPs) = lmT\,Pi) = Y,c^cit{n,m)\n,m). (3.2) 

n,m 

The probability for diffractive scattering is given by 

m^s) = Y,P^Pt[t{n,m)f = {t^)R,L. (3.3) 

n,m 

We have here identified |c„|^ = P„ for both R and L with the probability distribution 
for the dipole configurations inside the particles. Note that the sum Yin actually involves 
a sum over the dipole occupation number n as well as integrations over the transverse 
coordinates of each dipole, together with sums over their colour and spin configurations. 

Using the completeness of the states {\N,M)}, the expression in eq. ( ^ ) can be 
written in the following form: 

{i;s\i^s) = \{R,ms)\'+ \{N,L\i^s)\'+ ^ |(i?,M|^5)|' + 

+ \{N.Mm?. (3.4) 

Here the first term on the RHS corresponds to elastic scattering, where both R and L 
emerge intact from the collision. The second (third) piece gives the probability of the 
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excitation of R (L) into one of the states (M) with L (R) remaining intact. This 
corresponds to single diffractive excitation. Finahy, the last term takes into account the 
fact that both R and L may transform into excited states and M, which thus corresponds 
to double diffractive excitation. 

We note that the different terms in eq. ( |3.4D correspond to different averages of t{n, m)^. 
The sum of the single diffractive excitation and the elastic cross section can be calculated 
as follows 



N N 



I {N, L\i.s)\' = E E * = E E 'n* c^nt{n))i 

N n,m N n,n' 

Y.Pn{t{n))l = {{t)l)R (3.5) 



where we used the completeness of the states {A^}, 

^c^*c^ = 5„„,. (3.6) 



TV 



Each of the coefficients c„ above is to be evaluated at a certain rapidity Yq. The total 
rapidity interval between R and L is determined by the total cms energy ^/s of the process. 
For pp scattering Y is simply given by ln(s/Mp), while the situation is a bit more subtle 



in DIS. How we determine Y in DIS was discussed above in section 2.3.1. 

In Mueller's dipole model the scattering amplitude t(n, m) is given by the eikonal form 
1 — e~^, where F = ^ fij is defined in eqs. (^), (|2.3| ). The different contributions to 
the diffractive cross section in eq. ( |3.4D are then given by 



1 - e-^\ (3.7) 

/ R,L 

d-SD _//,_^-Fy\ _/l_e-^V (3.8) 



LI R \ / R,L 



SD _ / / ^ _ ^-fV \ _ /i _ ^^f"" ^ 



d^b 



Rl L \ I R,L 



i?,,L \ \ / LI R \\ I Rl L 



2 

^F ' 



+ ( 1 - e-M . (3.10) 

\ / R,L 

Here o^j-) (c^/j) is the cross section for the diffractive excitation of R (L). Similarly Of-i 
and Odd stand for the elastic and double diffractive cross sections respectively. Summing 
these four contributions we get the total diffractive cross section 



d^b 



{i'sHs) = Y.PnP^ [t{n,m)]' = /(l-e-A'\ . (3.11 

nm, \ ^ ^ ' R,L 
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Assume now that the state R is evolved up to Yq while L is evolved up to Y — Yq, with 
Y the total rapidity interval. The total and elastic cross sections given by 



atot{Y) = 2 1 d26^i_e-^(^)^ and a.iiY) = J cfh(l-e-^^^'>'j (3.12) 

are necessarily independent of Yq due to the requirement of frame independence. 

The diffractive cross section in eq. ( 3.11| ) is, however, not independent of Yq- This 



expression gives the probability for diffraction at a particular value of Yq, i.e. the chance 
that we find a rapidity gap around that particular Yq. If we calculate e.g. a^j^ at a specific 
Yq, we obtain the cross section where the diffr actively excited right-moving particle is 
confined within the rapidity range (0,yo)- This is approximately equivalent to a maximal 
diffractive mass given by M'j^ e^" . Taking the derivative with respect to Yq therefore 

gives the mass distribution da^j^/dYQ = dagj-)/dln{M^). 

For the total diffractive cross section in eq. ( p. Ill ) we thus expect, in the case of the 
symmetric pp collision, a maximum when Yq = Y/2, and a decrease when either Yq — > 
or Yq ^ Y. For asymmetric scattering as in DIS, we expect that it is easier to excite the 
photon. For a right-moving photon the diffractive cross section should therefore be smallest 
when Yq — > 0. This discussion is in accordance with the analysis by Hatta et al. [17]. In 
this reference the diffractive excitation of the proton is neglected, while in our formalism 
we can also take the proton excitation into account. 

The results obtained for pp collisions and DIS are presented in section ^. 

3.2 Importance of fluctuations 

The impact of fluctuations upon the small-x evolution has gathered considerable interest 
lately. As mentioned above, the various expressions for the cross sections in formulas 



eqs. (^?^)-( 3.10 ) are all obtained by taking different averages of the quantity (1 — e~^)^. 
The diffractive excitation is therefore completely determined by the fluctuations in the 
colliding systems and the interaction probabilities. 

3.2.1 Different sources 



There are several sources of fluctuations in the various expressions in eqs. (|3.7D -( ^TC ), 



related to variations in impact parameter, in the dipole cascades, and in initial the wave 
functions for the photons and protons. Many analyses include part of the fluctuations, 
assuming this to give the dominant contribution. Thus the dipole-saturation model by 
Golec-Biernat and Wiisthoff [7, 8] takes into account fluctuations in the photon wave func- 
tion and the emission of the first gluon in the photon cascade, while the model of Kowalski 
and Teaney (KT) [29] emphasizes the fluctuations in the impact parameter. Hatta et 
al. [17] includes the fluctuations in the photon cascade, but assumes that the fluctuations 
in the proton cascade can be neglected. As a result fits to data can give different results 
for the impact parameter profile, and different approaches can give different ratios for the 
elastic cross section and diffractive excitation. As an illustration we will here compare the 
fluctuations in our model with those in the Kowalski-Teaney model. 
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In the KT model the differential dipole-proton cross section in impact parameter space 
is given by the eikonal 

^ = 2(l-e-^('-''')/2), (3.13) 
where r is the dipole size. The opacity 0, is modeled by a factorized form 

^{r, b) = ^r^as{fi^)xg{x, y?)T{b), (3.14) 

where T{h) is the transverse profile function of the proton. In VL{r,b) there is also a 
dependence on W"^ oc 1/x which is omitted here. To determine the impact parameter 
profile it is assumed that the t dependence of the diffractive vector meson production cross 
section is given by an exponential, which in turn implies a Gaussian profile for T{h). The 
two unknown parameters of this Gaussian are then determined by a fit to diffractive J 
production data. 

For a virtual photon the only fluctuations are those in the dipole size, r, and impact 
parameter, h. The diffractive cross section is calculated as 

= j <fbj <frdz^l^^{r,z,Q^) [l - e-^M)/^^, (3.15) 

where ^^^{r, z, Q"^) = li^rir, z, Q^)^ _^ \^^{r, z, Q^)^. 

When comparing the two models it may seem natural to compare 0(r, 6)/2 with the 
average where {)d,p denotes the averaging over the dipole and proton cascades and 

the initial proton wavefunction. The corresponding total and diffractive cross sections 
would then read 



atot = 2J d'bj d'r dz ^^{r, z, Q') ^1 - e"^^ ^' '"^"'^ j , (3.16) 

= Jd'bJ d\dzi;,{r,z,Q') (l - e-^^^^'^^^^-^ . (3.17) 

This is, however, not necessarily correct. The opacity is in the KT model determined 
by a fit to data for the total cross section. What is directly determined is therefore 1—6"^/^, 
rather than Q itself. Thus a more direct analogy to our model would be the quantity 
(1 — e~^), which would give the same total cross section as eq. ( ^.12 ) and imply the 
following form for the diffractive cross section: 

a£}^ = I d'b I d2r.dzV7(r-,^,Q')^l-e-^(^'^)^' . (3.18) 

Note that this is not the same as the elastic contribution in eq. ( p.7[ ) since in that case also 
the photon wave function is included in the squared average (eq. ( |3.7[ ) is meaningless in 
DIS since the virtual photon cannot scatter elastically) . 

The expressions in eq. (3.17) and eq. ( p.lSj ) should be compared with the results in 
our model, obtained by integrating eq. ( 3.11 ) over b: 

^diff = Jd'bJd^rdzi:^{r,z,Q')(^(^l-e-^^'-'^^y^ . (3.19) 
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Figure 7: The ratio of tlie difFractive cross section to the total cross section for Mx < 32 GeV^. 
Our results obtained from cq. ( 3.19 ) (full lines) are compared to results obtained from both 
eq. ( 3.17 ), marked ((1)) in the left figure (dashed lines), and eq. ( 3.1g ), marked ((2)) in the right 
figure (dotted lines). The total cross section is calculated according to eq. (B.12). 



Such a comparison is interesting as a way to gauge the role played by the fluctuations. 

The results obtained from eq. ( |3.17| ) and eq. ( |3.18| ) are shown in figure 0, together with 
the results obtained from eq. ( 3.19 ) which includes all fluctuations. We immediately notice 
the very large effects in our model from the fluctuations in the cascades and the proton 
wave function. The diffractive cross section calculated from eq. (|3.17] ) is seen to be a factor 
2-3 below the result obtained from eq. ( p. 19 ), while the result from eq. ( |3.18| ) is around a 
factor 4 lower than eq. (|3.19| ). 

We conclude that in our model a large fraction of the fluctuations determining the 
diffractive cross sections is caused by the dipole cascade evolutions. In order to obtain 
a similar result in the KT model it is therefore necessary to have larger fluctuations due 
to the impact parameter dependence, which means an impact parameter profile which is 
more narrow, i.e. more black and white compared to Mueller's dipole cascade model, where 
the average scattering can be "grey" overall, since the fluctuations in the cascades means 
that some events are almost black while other are almost white. This can clearly be seen 
in figure |8| where we compare the impact parameter profile from the KT model to that 
obtained from our model. (A similar effect, although less pronounced, is observed in the 
profile for pp scattering presented in ref. [15]. Also here the profile obtained in our model 
has a somewhat higher tail for large impact parameters than the Gaussian fit to Tevatron 
data by Sapeta and Golec-Biernat [30].) 



3.2.2 Wave functions 

Photon 

The photon wave functions in eq. (^]^) for longitudinal and transverse photons are fully 
determined by perturbation theory. We note that these wave functions are not normalized, 
even for real photons. This is, however, not in contradiction with the assumption that the 



-17- 




Figure 8: The impact parameter profile of dipole-proton collisions for an initial dipole size 
r — 2GeV~^ ~ 0.4fm at two different energies, y = Ins = 9 and 11. Om- results (solid and dashed 
lines) are compared to those from the Kowalski-Teaney (KT) [29] model (dot-dashed and dotted 
lines). 



states {|A^)} are normalized. The generic photon state can be written 

oo 

|7(yo)) = cliYo) |7d) + Y.'^n{Yo) In), (3.20) 

n=l 

where [7^) is the component of the photon coupling directly to the quarks. While the 
state I7) is normalized, the two separate components above are not. Rescattering of the 
component [7^) can be neglected, as it is proportional to aem- It is therefore sufficient to 
only keep the contribution from \ jh) = ^n^n\n), which is not a normalized state. 

In DIS it is not meaningful to consider the elastic j*p scattering, as the virtual photon 
can never be detected as a real particle. The closest analogies to elastic scattering are given 
by Deeply Virtual Compton Scattering (DVCS) and the exclusive reactions 7*p Vp, with 
V a vector meson. We will return to these processes in a future publication. 

Proton 

The wave function for the proton is much less well defined. The expressions in eqs. ( ^.^ ) 
or (^^) ought to be interpreted as probability distributions rather than quantum mechan- 
ical wave functions, which can be used to determine the interference effects present in 
eqs. (|3^)-(|3lO|). The fluctuat ions in the wave functions influence the terms in eqs. 



(3.10) in which the average of 1 — e is taken before the square {i.e. events containing an 



elastic proton). 

As discussed in section 2.3.2 a wave function where three particles can simultaneously 
be in a single point is not realistic. A Gaussian distribution in impact parameter is usu- 
ally motivated by the exponential dependence on t for the elastic cross section. However, 
the constant t-slope is (except possibly for the highest Tevatron energy) only valid for 
\t\ < 0.15 GeV^, corresponding to 6 > 2/\/0.15 GeV~^ Ifm, and therefore a suppres- 
sion for small r is still compatible with this constraint. Simulations with such a wave 
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function reduces the fluctuations and increases the cross sections for elastic scattering and 
single diffractive excitation, while leaving the total and the total diffractive cross sections 
unchanged, provided the average (r-^) is kept the same. 

There is, however, also a more fundamental problem with the proton wave function. 
In the Good and Walker formalism the hadronic states {|iV)} form a complete set. This 
implies that before the cascade has started, there must be other hadron states with wave- 
functions orthogonal to the proton wavefunction. This calls for a detailed dynamical scheme 
describing the relevant degrees of freedom for the hadronic states. With the approxima- 
tion expYb these orthogonal states also have the same mass as the proton, which 
increases the problem further. 

Lacking a real quantum-mechanical description of the proton wave function, we can still 
get an upper limit for elastic scattering and single diffraction by removing the contribution 
from the initial wave function fluctuations. This is obtained if we integrate over the initial 
wave functions in eqs. (|3.7|) - (|3.10| ) after taking the squares. Note that the average over 
different evolutions is still taken before the square, and therefore the fluctuations in the 
cascade evolution and the impact parameter dependence are still included. Note also that 
this does not affect the result for the total cross section in eq. ( p.l2D or the total diffractive 
cross section in eq. ( |3.11| ) (which also includes the elastic cross section). 

3.2.3 Non-leading eflfects 

It was early pointed out by Mueller and Salam [31] that there are extremely large fluctu- 
ations in the leading order cascade evolution. Expanding the exponential in (1 — e~^) we 
have 

l-e-^) = Y.^-^{F'). (3.21) 

' k=l 

Here {F^) could be interpreted as a contribution from the exchange of k pomerons, but 
such an interpretation may be difficult as it was numerically demonstrated^ by Salam [32] 
that {F^) ~ (A;!)^, and therefore this series is strongly divergent. The reason for this is the 
existence of rare events with a large number of dipoles and large values of F, which make 
{F^) blow up for large k. On the other hand, these rare events do not contribute much 
to (1 — e~^), since this expression saturates for large F. Due to these large fluctuations 
it is therefore possible that one can observe events in which there are large saturation 
effects, even though the average scattering is still weak. Although such rare events are 
less important for the total cross section, they are very essential for the diffractive cross 
section. 



In leading order the dipole splitting in eq. (2T) diverges for small dipoles, and therefore 
the number of small dipoles depends strongly on the necessary cutoff. As pointed out in [32] 
a very essential source of the large fluctuations is also the occasional creation of a very large 
dipole. Such a large dipole has a large probability to split, and the most likely scenario is 



^The numerical result was anticipated by Mueller [11] who performed analytical calculations on a simple 
toy model in which transverse coordinates are neglected. 
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that it splits into one very small dipole and one dipole which is almost equally large. The 
process is then iterated and the result is a "jet" of many small dipoles. 

We note, however, that the very large fluctuations observed by Mueller and Salam 
are strongly reduced by non-leading effects. As demonstrated in [14] energy-momentum 
conservation has a very strong influence. The production of small dipoles is suppressed 
by the conservation of the positive light-cone momentum while the large dipoles are 
suppressed by the conservation of p- . Another non- leading effect comes from the running 
coupling Og. As discussed in section the relevant scale in the dipole splitting is de- 
termined by min(r, ri,r2) where r is the mother dipole which splits into ri and r2. This 
suppresses not only the production of very small dipoles, but also fluctuations where a 
small dipole splits producing two very large dipoles. We thus conclude that both effects 
contribute to a suppression of very small or very large dipoles, and therefore also of the 
abovementioned "jets" radiated from an occasional large dipole. As a result we find in 
our calculations that {F ) grows like kl rather than (A;!)^. The ratios {F^)/{k ■ {F^~'^)) for 
dipole-dipole scattering at y = 10 are approximately equal to 1.2 for all k between 5 and 9 
(larger /c- values need very high statistics, and are therefore difficult to simulate). However, 
although the fluctuations in the cascade evolution are strongly tamed by non-leading ef- 
fects, they are still very important, and have a large effect on the diffractive cross sections, 



as seen in section p.2.l| and in the results presented in the next section. 



4. Results on diffraction and elastic scattering 
4.1 Diffraction in pp collisions 

In figure |^ we show the ratio of the total diffractive, single diffractive and elastic cross 
sections to the total cross section at 1.8 TeV as a function of the Lorentz frame used. 
Figure |^ shows the result obtained with the initial proton wavefunction in eq. ( |2.8| ) , while 
figure shows the upper limit obtained by integrating over the initial proton wavefunctions 
after taking the square in eqs. (3.7)-(|3^). In these figures we also show results from the 



CDF collaboration (not including double diffraction). 

The results in figure ^a do not agree well with the data. The elastic cross section is 
too low. The single diffractive cross section (T^^, does not go to zero when Yq — > 0. In this 
limit M|- nip and there should be no phase space for diffractive excitation. The double 
diffractive excitation is about 9%, which is unrealistically high, and this cross section should 
also approach zero in the limits Yq^OoiYq^Y. All these features illustrate the problems 
with our initial proton wave function, discussed above. This formalism presumes that there 
are other hadronic states, which have wavefunctions orthogonal to the proton wavefunction, 
and which have approximately the same mass. The limiting results in fig obtained when 
the fluctuations in the initial proton wave function are neglected, do not have this problem. 
We also see that these results agree quite well with the experimental estimates, supporting 
the assumption that the initial wavefunction fluctuations have a small effect compared to 
the fluctuations in the cascade evolution and the impact parameter dependence. 

For the elastic cross section we also note (as was claimed above) that it is approximately 
frame- independent, as it should be. The upper limit shown in figure lb is around 22%, 
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Figure 9: The ratio between tlie total difFractive and the total cross section (solid line) together 
with the contribution from elastic (dashed), single-right (dash-dotted) and double difFractive (dot- 
ted) cross sections at f .8 TeV, obtained including (a) and excluding (b) fluctuations in the initial 
proton wavefunction. In both cases the lower error band is an estimate from CDF data on single 
diffraction [33], the middle band is the CDF elastic cross section [34] and the upper is a sum of 
the two contributions. (Thus the contribution from double diffraction is not included in the CDF 
result.) 



which agrees well with the value (22.02 it 0.78)% from E811 [35], while the value from CDF, 
(24.6 ± 0.4)% [34] is a bit larger. 

The single diffr active cross section a^j^ (^sd) increasing (decreasing) when Yq is 
increased, and in the model we make the identification of the diffractive masses: 

M^(^) = GeV2, M^^^^ = e^-^« GeV^, (4.1) 

where for 1.8 TeV we have Y = ln(s/lGeV2) w 15. In fi gure ^ we also show results obtained 
from the CDF parameterization of single diffractive excitation [33] : 

dcrf^ _ 1 D f s 

d\nM^^^^ ~^{hQ + 0.51n(s/M^^^^)) VmJ^^ 

with D = 2.54 ± 0.43 mb, 6o = 4.2 ± 0.5, and e = 0.103 ± 0.017. Similar results are 
also presented by the E710 collaboration [36], although with somewhat larger errors. Our 
results in figure ^b agree quite well with the data, even if they are a little high for the 
largest excited masses. Besides not going to when Iq ~* Oi the result in figure ^ also 
has a much too slow variation with Yq, meaning a too low value for dasD / dhiM'j^. 

For double diffraction our result for Yq = 7.5, which corresponds to a central gap, is 2.0 
mb. Experimental data exist for 900 GeV from the UA5 collaboration at the CERN SppS 
collider [37]. Our result at this energy is 1.8 mb, which is consistent with the experimental 
result 4.0 it 2.2 mb. Our results can also be compared to the model of Goulianos [38], 
who argues that aoD should decrease with energy, due to saturation effects, from around 
1.6 mb at 900 GeV to around 1.3 mb at 1.8 TeV. We also note that our result is consistent 
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Figure 10: The ratio of the elastic to the total cross section in collisions as a function of ^fs for 
various values of r,„aa:(GeV~^) and AQC£)(GeV). Lower curves are obtained including fluctuations 
in the initial proton wavefunction, while the upper curves excludes these fluctuations. The data 
points arc from CDF [34] and E811 [35]. 



with a factorized dependence on the two masses, as expected from Regge formalism: 

where /(M^) denotes the distribution for single diffraction in eq. ( [4. 2D . (Note that double 
diffraction is not included in the CDF data in figure |9|.) 

As seen in figure § the same total cross section can be obtained with different sets 
of values for Kqcb and Vmax- In figure ^ we see that varying these parameters, keeping 
Gtot constant, does modify the elastic cross section somewhat. However, our upper limit is 
still close to the data, leaving little room for a contribution from the initial proton wave 
function. 

Finally we show in figure |ll] the energy dependence of the total, the total diffractive 
(including elastic scattering), and the elastic cross sections, including our predictions for the 
LHC. (The values presented correspond to those in fig ^3, i.e. to the more realistic results 
obtained neglecting the fluctuations in the initial proton wave function.) The diffractive 
cross section is calculated in the cms with Yq = Y jl, which demands a central gap and 

predict for the LHC are 108, 38, and 26 mb respectively. 



implies that the diffractive excitations are limited by M\ < y/s ■ IGeV. The values we 



4.2 Diffraction at HERA 

Diffractive excitation has been measured at HERA by the ZEUS [39] and HI [40] collabo- 
rations with two different methods. One is based on an observed rapidity gap. The ZEUS 
data obtained with this method [39] give the cross section integrated over all diffractively 
excited protons with mass < 2.3 GeV. Assuming that the contribution from events 

where the proton is excited beyond this limit is small, and can be neglected, the result 
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Figure 11: The total (full line), total diffractive (including elastic) (dashed line) and elastic cross 
sections (dotted line) as a function of collision energy, -^5, in pp collisions. Here the diffractive cross 
section is evaluated at Yq/Y = 0.5 i.e. with a central gap in the cms. 



of this method for m]^-* < M^^'^^^ corresponds to our model calculations for (J^i// at 
Yo = ln(M 



r(7) 



2(7) ^ 
X,max ) 



InM 



2(7) 
X,max 



d\nM 



2(7) 
X 



da. 



diff 



dlnM 



2(7) 
X 



-Sr^^o = ln(MS.J). (4.4) 



The results are shown in figure 12 and we see a very good agreement with data, although 
there is a tendency for our cross sections to decrease a bit too slowly with Q^. 

The cross section for single diffractive excitation of the photon can also be calculated 



in our model, and in figure 13 we present the ratio wrt. the total diffractive cross section 
as a function of W for different and Mx ■ In [39] the ZEUS collaboration estimated this 
ratio to be 0.70 it 0.03, by comparing a parameterization'^ of their diffractive data to results 
from their leading proton spectrometer. This result is obtained using the assumption that 
the ratio is independent of W, and Mx- Comparing with figure 12 we find that our 
result is consistent with the ZEUS number, but that we predict that the ratio actually does 
have a small dependence on Mx and Q^. 



5. Conclusions and Outlook 

We have in earlier papers presented an extension of Mueller's dipole cascade model, which 
includes subleading effects from energy conservation and running coupling as well as colour 
suppressed effects from pomeron loops. The model is also implemented in a MC simulation 
program, which simplifies the comparison between theoretical ideas and experimental data, 
and allows more detailed studies of important non-leading effects. Calculations of total 

^Using a modified version of tfie model in [41]. 
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Figure 12: The ratio of the total difFractive cross section to the total cross section as a function 
of W, for Mx < 4 (a), 8 (b), 16 (c) and 32 GeV (d). Our results are compared to ZEUS data [39] 
for = 4 (full lines and + points), 14 (dashed lines and x points) and 55 GeV^ (dotted lines and 
open squares). 

cross sections agree very well with experimental results for pp collisions and deep inelastic 
electron scattering. 

To gain further insight into small-x evolution and saturation, we have in this paper first 
presented a number of improvements of the model, in particular related to the confinement 
mechanism, and thereafter applied the model to elastic scattering and difFractive excitation, 
where we specifically study the effects of different sources for fluctuations. 

A consistent treatment of confinement effects is achieved by replacing the infinite range 
Coulomb potential in the dipole splitting and in the dipole-dipole scattering with a screened 
Yukawa potential. By equating the screening length, rmax) with the size of the proton 
entering into its wavefunction, we were able to get a good, boost-invariant description of 
the pp and 7*p total cross sections for a wide range of energies, using basically only two 
parameters, rmax and Aqcd- This new treatment of confinement has effects on the boost 
invar iance of the model, further improving the earlier, almost frame independent, results. 

Our treatment of diffraction is based on the formalism of Good- Walker and Miettinen- 
Pumplin. The cross sections for elastic scattering and diffractive excitation are here deter- 
mined by the fiuctuations in the interaction probability between different events. Contrary 
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Figure 13: The ratio of the single difFractive cross section to the total diffractive cross section as 
a function of W, for Mx < 4 (a), 8 (b), 16 (c) and 32 GeV (d) for = 4 (f^n lines), 14 (dashed 
lines) and 55 GeV^ (dotted lines). 



to other calculations, we can in our model easily consider all different sources for such fluc- 
tuations; those stemming from the dipole cascade evolution, the dipole-dipole scattering, 
from the impact parameter dependence, and from the initial photon and proton wavefunc- 
tions. We find that all of these sources give important contributions, apart from the initial 
proton wavefunction, and together they give a very good description of data on elastic and 
single- and doubly-diffractive scattering in both 7*p and pp collisions. We must, however, 
admit that wc do not have a realistic quantum-mechanical description of the proton state 
in terms of dynamical variables. Here data are best reproduced if the contribution from 
the fluctuations in the initial proton state are small compared to the other contributions. 

In a future publication we will use our model to study the quasi-elastic reactions 
7*p Vp and deeply inelastic Compton scattering (DVCS). In the future we also want 
to develop the model further to be able to describe exclusive multi-particle final states. 
This needs, however, a recipe for how to handle the virtual dipoles, those which do not 
participate in the collision and therefore cannot come on shell and give final state hadrons. 
To describe particle production in the proton fragmentation regions would also need a 
much improved description of the initial proton state. 
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